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Abstract. We present a detailed chemical analysis on 
the basis of high-resolution, high signal-to-noise optical 
spectra of two post-AGB objects IRAS 22223+4327 and 
IRAS 04296+3429. Both display the unidentified 21 
feature in their IR-spectra. The analysis is performed 
using HR 1865(F0Ib) and HR 1017(F5Iab) as reference 
stars. 

The spectroscopic indicators provide accurate atmo- 
spheric parameters of T e ff=6500 K, log</=1.0 and 6=5-5 
kms" 1 for IRAS 22223+4327 and T off =7000 K, log #=1.0 
and 6=4.0 kms" 1 for IRAS 04296+3429. Our high- 
resolution data are inconsistent with the significantly 
lower temperatures deduced from spectral-type determi- 
nations in the literature based on low resolution spectra 
and highlight the need of high-resolution spectroscopy for 
the determination of accurate fundamental parameters of 
chemically peculiar supergiants. 

Both photospheres are found to be metal-deficient with 
[Fe/H]=— 0.4 and —0.7 respectively. C and N are found 
to be overabundant. Useful O-lines were only detected 
in the brighter IRAS 22223+4327 and the O abundance 
is found to follow the Fe deficiency : the C/O photo- 
spheric abundance is about 1.3, but due to the lack of 
oxygen lines it is difficult to determine accurately. This 
corroborates the fact that the carriers of the 21 /zm fea- 
ture are formed in a carbon rich circumstellar chem- 
istry. Moreover, these IRAS-stars have large overabun- 
dances of s-process-elements. The mean abundance of 
all the measured s-process-elemcnts is [s/Fc]=+1.0 for 
IRAS 22223+4327 and +1.4 for IRAS 04296+3429. The 
distribution of the s-process elements can best be de- 
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scribed as due to a distribution of neutron exposures with 
a low mean neutron exposure of To = 0.2 mbarn -1 . The 
21 /j,m stars form an interesting sub-group in the total 
post-AGB sample of stars, not only for their IR charac- 
teristics, but also in a broader context of stellar (chem- 
ical) evolution theory. They show, in contrast to other 
post-AGB stars, that the 3rd dredge-up has been efficient 
during their AGB evolution. The mean neutron exposure 
is lower than expected for their metallicity. 

The spectroscopic parameters found for the massive 
spectral analogues (T c ff=7500 K, logg=2.0 and 6=3. 
kms" 1 for HR 1865 and T cff =6500 K, logg=2.0 and 
6=3.5 kms -1 for HR 1017) agree well with the values 
found in the literature. Both stars display an enrichment 
in Na, which is commonly observed in massive supergiants 
but theoretically not well understood. 

Key words: stars: individual: HR 1017, HR 1865, 
IRAS 04296+3429, IRAS 22223+4327 - stars: post-AGB 
- stars: abundances 



1. Introduction 

The late AGB and post-AGB phases are among the least 
understood in current low- and intermediate-mass stellar 
evolution theory, mainly because of the coarse understand- 
ing of the mass-loss mechanism(s) and internal convection. 
The theoretical understanding of the total chemical evo- 
lution and certainly the chemical enrichment due to the 
3rd dredge-up phenomena are therefore still poorly under- 
stood (e.g. Lattanzio et al. 1996). 

One of the tracers of the chemical evolution of the star 
is the chemical content of the dusty circumstellar envelope 
which is formed by the high mass-loss episode (s) during 
the late AGB phase. The infrared spectral region is often 
used to infer the envelope composition since not only the 
bulk of radiation is emitted in the IR but also because 
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the IR-spectra are characterized by the chemistry. The 
main signature of O-rich dusty environments are silicate 
features at 9.7 /im and 18 /an while a C-rich chemistry is 
characterized by the presence of the SiC emission at 10- 
12 /im and/or the presence of 3.3, 6.2, 7.7 and 11.3 /im 
features which are usually attributed to the polycyclic aro- 
matic hydrocarbons (PAHs) . The spectrographs on board 
of ESA's Infrared Space Observatory (ISO) are revealing 
now a wealth of detail concerning these 2 chemical types 
of circumstellar environments around evolved stars (see 
e.g. the ISO dedicated November 1996 issue of A&A). 

Recently the detection of a strong broad emission fea- 
ture at 21 /mi accompanied by unusually strong 3.4-3.5 
and 6-9 /im emission features in several post-AGB stars 
has been reported (Kwok et al. 1989; Hrivnak & Kwok 
1991a,b; Kwok et al. 1995; Geballe et al. 1992; Henning 
et al. 1996). While the carrier of the broad feature is not 
yet identified (e.g. Begemann et al. 1996 and references 
therein), there are several observational indications that it 
is only present in a C-rich circumstellar chemistry: the 21 
/im stars usually also display strong PAH features at lower 
wavelengths; Hrivnak (1995) studied low resolution spec- 
tra and revealed the presence of circumstellar C2 and C3 
and also Bakker et al. (1995,1996,1997) detected the op- 
tical circumstellar carbon molecules around 21 /im stars; 
Omont et al. (1993) deduced a C-rich environment from 
the large HCN/CO millimetre-line ratio. Till now the 21 
/im feature is only observed in post-AGB stars and young 
PN's and it is not clear whether the carrier is only pro- 
duced during the transition from the AGB to the PN phase 
or already produced on the AGB but not excited. 

Since the discovery of Justannont ct al. (1996) that 
also HD 187885 may display a small 21 /xm feature, all 
known the post-AGB stars with photospheres showing the 
yields of an efficient 3rd dredge-up (Van Winckel 1997 
and references therein), also display the 21 /im feature in 
their IR spectrum. The study of the stellar photospheres 
of the 21 /im stars is therefore not only interesting for the 
study of the carrier and excitation mechanism of the 21 
/im feature itself, but also in a broader context of stellar 
(chemical) evolution of low and intermediate mass stars. 

In this paper, we report on a detailed chemical study 
of the stellar photospheres of two 21 /im post-AGB stars: 
IRAS 22223+4327 and IRAS 04296+3429. The two ob- 
jects are optically faint (V = 9.6 and 14.2 respectively) 
with a total optical flux which is smaller than the total 
infrared flux. Their spectral energy distribution is double- 
peaked, with visible and near-infrared components due 
to the reddened photosphere and a far-infrared compo- 
nent caused by the detached dust shell (Hrivnak & Kwok 
1991a). The huge IR excess may indicate a relatively high 
initial mass. In the IRAS colour-colour diagram, the two 
stars are located in the region of post-AGB stars in be- 
tween the late AGB stars and young PN's (Omont et al. 
1993). We used two bright massive supergiants HR 1017 
(a Per) and HR 1865 (a Lep) of similar spectral type 



(F5Ib and FOIb respectively) as reference objects (see Ta- 
ble 0). 



2. Observations and data reductions 

The chemical analysis of the programme stars is based 
on an accurate study of high-resolution optical spectra 
(A/<5A=50000, corresponding to a projected slit of 1.1"). 
The spectra have been obtained by one of us (EJB) with 
the Utrecht Echclle Spectrograph, fed by the 4.2m William 
Herschel Telescope at La Palma, Spain. The echelle spec- 
trograph is equipped with a cross-dispersor so that the 
different orders are projected next to one another on the 
CCD. In the red part of the spectrum the consecutive or- 
ders are not overlapping. The wavelength coverage of the 
spectra is indicated in Table |l} The data were reduced us- 
ing the echelle package of IRAF V2.9 in a standard way 
with optimal extraction of the orders. 

3. Analysis 

3.1. Atomic Data 

For the line identification wc mainly used the line lists 
by Thevenin (1989, 1990), which are based on the lines 
identified in the solar spectrum by Moore et al. (1966). 
These lists are completed with the line lists of Van Winckel 
(1995) and Bakker (1995). We first performed a complete 
identification of the high S/N spectra (S/N > 250) of 
HR 1865 and HR 1017 which we used to identify and de- 
blend the lower quality spectra of the IRAS sources. 

In recent years considerable efforts have been made in 
calculating improved model stellar atmospheres and os- 
cillator strengths (\oggf) for atomic transitions. Unfor- 
tunately, high precision gf-values are still unavailable for 
many elements. We have taken the excitation potential 
and the oscillator strengths mainly from the inverted solar 
abundance analysis of Thevenin (1989, 1990). The atomic 
data for O I, Mg I, Mg II, Si I, Si II, Ti I, Ti II, Cr I, Cr II, 
Fe II, Sr I, Sr II, Y I, Y II, Zr I and Zr II are obtained 



from the automatic data bank of Kurucz (http://cfa- 



www.h~arvard.edu/amp/data/kur23/sckur.htmj ) and the 
mailserver of the 'Vienna Atomic Line Data Base'. For 
the Fe I oscillator strengths we used the critical compila- 
tion of Lambert et al. (1996). 

3.2. Determination of the atmospheric parameters 

We used the CDROM grid of LTE model atmospheres 
of Kurucz (1993) in combination with his abundance cal- 
culation programme WIDTH9. A model atmosphere is 
uniquely determined by the metallicity, the effective tem- 
perature (T e ff), the gravity (log g) and the microturbulent 
velocity (£t). We used input models with a solar metal- 
licity Z for HR 1865 and HR 1017 and Z=-0.5 for the 
IRAS-sources with £t=2 kms -1 . 
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The determination of these atmospheric parameters 
were solely based on our high-resolution spectra. Quan- 
titative photometric analysis is hampered by the uncer- 
tainty on the amount of circumstellar reddening, the pos- 
sible anomalous circumstellar reddening law and by the 
inaccurate calibration of the photometric systems for su- 
pergiants in general. The photometry can therefore only 
be used as a first guess for the finer spectroscopic analysis. 

3.2.1. The effective temperature 

A spectroscopic estimate for T e f f is found by forcing the 
abundances of an ion to be independent on the lower exci- 
tation potential of the transitions. For F-type supergiants, 
only Fe I lines are suitable for such an analysis since this 
method is only usable for ions with a large number of 
weak lines having a large spread in excitation potential. 
Moreover, for supergiants in the temperature range of our 
programme stars, non-LTE effects on the Fe I lines are 
small (Venn 1995b). 

During our research it became clear that even a few un- 
precise oscillator strength values can limit the precision of 
this method considerably. When using the Fe I oscillator 
strengths of Thevenin we encountered the same problem 
for every star : the effective temperature determined on 
the base of the Fe I lines was some 500K higher than the 
temperature based on the Fe II lines. Moreover, the tem- 
perature estimates of the reference stars were higher than 
the temperatures found in the literature. 

First we situated this discrepancy in the departure of 
local thermodynamic equilibrium (non-LTE effects). The 
LTE underabundance of iron due to overionisation of Fe I 
in F supergiants amounts from ~0.03 dex (F8 supergiants) 
to 0.2 dex (F0 supergiants) (Venn 1995b; Boyarchuk & 
Lyubimkov 1983; Boyarchuk et al. 1985, 1988) A shift of 
0.2 dex in abundance corresponds to a temperature shift 
of ^200 K. It turned out, however, that for every star 
AT eff (=T eS (Fe I)-T cff (Fe II)) was too large (e.g. 700 K 
for HR 1865 and 400 K for HR 1017) to be caused solely 
by non-LTE effects. Note that by limiting our abundance 
analysis to lines with an equivalent width smaller than 150 
mA, we focus on lines formed in deeper layers where the 
non-LTE effects are known to be smaller. 

We have then limited the Fe I lines to only those listed 
in the critical compilation of log (//-values by Lambert et 
al. (1996). They re-discussed the Fe I gf-values by com- 
paring the gf-values of several authors (see Lambert for 
references). By using these oscillator strengths, all prob- 
lems concerning discrepancy of T e g, logg or £ t (non-LTE 
effects) disappear! It is extremely important to limit the 
used Fe I lines to transitions with a precise oscillator 
strength. For Fe II we don't have any lines in common 
with Lambert. Therefore we have compared the \oggf- 
values of Thevenin with these of Lambert: the \oggf- 
values of Thevenin were systematically 0.19±0.08 smaller 
than these of Lambert, due to the fact that Lambert used 



7.51 for the solar Fe abundance and Thevenin 7.67. We 
have thus increased the log gf- values of Fe II of Thevenin 
with 0.19. The solar abundances of the other chemical el- 
ements are from Grevesse (1989). For HR 1865, HR 1017 
and IRAS 22223+4327 we derived a temperature of 7500 
K, 6500 K and 6500 K respectively (e.g. see Fig. Q). 
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Fig. 1. The excitation potential-abundance diagram and the 
equivalent width-abundance diagram for IRAS 22223+4327. 

For the faint object IRAS 04296+3429 it turned out to 
be impossible to determine the effective temperature with 
this method due to the lack of good quality Fe I lines. A 
rough estimate is found when comparing the spectra of the 
IRAS-stars with those of the two reference stars. Figure 
^| shows the spectra of the four programme stars in the 
wavelength coverage between 6110 A and 6160 A. This is 
an interesting spectral interval since besides the O I triplet 
and Ba II line, also Fe I and Fe II lines are present. The 
ratio of the Fe II line-strengths versus the Fe I lines gives 
an indicator of the effective temperature for these four 
stars of similar gravity. 

The resemblance between the spectra of HR 1017 and 
IRAS 22223+4327 is striking. This is confirmed by our 
spectroscopic temperature assessment of T e ff=6500 K for 
HR 1017 and IRAS 22223+4327. For IRAS 04296+3429 
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the ratio of the Fe II to Fe I lines is higher than the ratio 
of HR 1017 and lower than this of HR 1865. We therefore 
estimate the temperature of IRAS 04296+3429 to be ap- 
proximately 7000 K. To test the consistency of these T e g, 
we have checked the excitation potential-abundance dia- 
gram of IRAS 04296+3429. The slope is small and pos- 
itive, but it indicates that the temperature departure is 
still less than 300 K, within the bounds of the estimated 
accuracy. 
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Fig. 2. Spectra of HR 1865, HR 1017, IRAS 22223+4327 and 
IRAS 04296+3429. 



Hrivnak (1995) has determined the spectral type of 
the IRAS sources on the basis of low resolution spec- 
tra and concluded that both IRAS 04296+3429 and 
IRAS 22223+4327 are GOIa supergiants. This spectral 
classification was used by Kwok et al. (1995) and Bakker 
et al. (1996) to estimate the effective temperature of 5000- 
5500 K for both stars. This low temperature is, however, 
inconsistent with our high-resolution spectra! Not only be- 
cause in the lower excitation potential- abundance plot, a 
significant upward trend shows up, but also because neu- 
tral lines of heavy elements should appear in the spectra 
at those temperatures. The strong Y I line at 6435 A for 
instance, would have an equivalent width of 20 mA for the 
abundance computed from the Y II lines at a temperature 
of 5500 K and log(g) = 1.0. 

We want to stress the fact that for reddened, chem- 
ically peculiar supergiants, high-resolution data are def- 
initely needed for accurate fundamental parameter de- 
terminations. The line regions, used in the low resolu- 
tion spectra to infer the spectral type, are influenced 



by strong resonance lines of s-process elements like Ba 
and Y. With strong s-process enhancements, like observed 
in IRAS 05341+0852 by Reddy et al. (1997), the opti- 
cal spectrum is even completely dominated by lines of s- 
process isotopes. In our opinion, the chemical peculiarity 
and certainly the strong enhancements of s-process iso- 
topes, make standard spectral classification difficult. 

3.2.2. Gravity and the microturbulent velocity 

We derived the model gravity by implying that different 
ions of the same element yield the same abundances to 
within 0.1 dex. For F supergiants, the only element use- 
ful for this purpose is again Fe by the lack of alternatives 
with enough lines of both ionisation stages. The observed 
hydrogen Balmer lines are affected by emission so we did 
not use them to constrain the gravity. A change in the 
logg value by 0.5 makes the abundances derived for the 
two ionisation stages differ by in between 0.11 and 0.26 
dex. The total uncertainty on the gravity of these objects 
is, however, more uncertain than the ± 0.3 this method 
suggests, since it does not take into account other un- 
certainties on the abundances which are more difficult to 
quantify, like non-LTE effects, systematic log gf errors etc. 
Note that the main conclusions of this work (see discus- 
sion) are based on abundance ratios, which are much less 
influenced by uncertainties in effective temperature and 
gravity than absolute values. 

The microturbulent velocity for each star has been 
determined by forcing the Fe I abundances to be inde- 
pendent on the equivalent width (Wa) (see Fig. |l| for 
IRAS 22223+4327). 

IRAS 04296+3429 has not enough lines available with 
a good range in equivalent width for a reliable estimate. 
We assumed £t=4.0 kms" 1 . Nevertheless, the value of mi- 
croturbulence is not critical in this analysis since only 
weak lines (Wa <150 mA) are included in the final av- 
erage abundances. The use of weak lines also means that 
radiative damping terms will not be important. 

A synopsis of the atmospheric parameters for the four 
stars is listed in Table |l|. 

3.3. Error analysis 
3.3.1. Internal errors 

The standard deviation a on the abundance of the ele- 
ment for which we observed more than five lines is a good 
indicator of the consistency of the chemical analysis (e.g. 
see Table |^). A typical value for a is between 0.10 and 
0.25 dex. 

For a good model in ionisation and excitation equi- 
librium, a is mainly determined by non-systematic errors 
on the equivalent width and especially on the oscillator 
strength. The equivalent widths were measured by fitting 
a Gaussian curve to the absorption lines. 
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Table 1. The atmosphere models for HR 1865, HR 1017, 
IRAS 22223+4327 and IRAS 04296+3429. 





m(v) 


T <sJf 




it 


Obs. region 




magn. 


K 


cm S -2 


kms" 1 


A 


IRAS 
04296+3429 


14.2 


7000 


1.0 


4.0 


5560-10045 


IRAS 
22223+4327 


9.7 


6500 


1.0 


5.5 


4430-10045 


HR 1017 


1.7 


6500 


2.0 


3.5 


4690-7340 


HR 1865 


2.6 


7500 


2.0 


3.5 


4690-7340 



To check our log (//-values for systematic errors, we 
compared our LTE abundance analysis of HR 1865 with 
a similar study of Venn (1995b). Venn mainly used the 
oscillator strengths of Fiihr et al. (1988); Wiese & Mar- 
tin (1980); Fiihr et al. (1981); Wiese & Fuhr (1975) and 
O' Brian et al. (1991). No systematic differences occur, 
so we continued to use the log ^/-values of Thevenin for 
most atomic species. Since we did not account for the hy- 
perfine broadening of transitions of atoms with an odd 
atomic number, the abundances of these ions are slightly 
overestimated. 

Other contributors to the scatter could be differential 
non-LTE effects and non-detected blends. The blends of 
the IRAS stars, which have broader lines, are detected 
by comparing the spectra with the narrow lined reference 
stars. 



3.3.2. Inaccuracies on the model parameters 

The model parameters (T e ff, log g and £ t ) are not inde- 
pendent: a change in one parameter generally induces a 
shift in another for the spectroscopic requirements (ioni- 
sation balance, independence of the abundance of an ion 
versus the excitation potential and equivalent width) to be 
fulfilled. A typical shift of 0.5 dex in the gravity induces a 
temperature shift of 300-400 K, so that the ionisation bal- 
ance still would be maintained. We refer to Table 3 of Van 
Winckel (1997) for a quantitative estimation of the influ- 
ence of the different uncertainties on the abundance de- 
termination. The temperature uncertainty has by far the 
biggest influence on the abundance accuracy, especially 
for ions with the smallest occupation level since they are 
most influenced by uncertainties on the ionisation balance 
induced by errors on the temperature and gravity. 



4. Results of the chemical analyses 

In the Tables |-| we give the results of the abundance 
analysis for each programme star. For every ion we list 
the number of lines used, the mean equivalent width, the 
absolute abundance, the abundance ratio relative to the 
solar value and the internal scatter, if more than one line 
is used. A complete line list with the detailed atomic data 
can be obtained upon request. For the solar iron abun- 
dance we used the meteoric iron abundance of 7.51. 



Table 2. Chemical analysis of HR 1865. In the last column we 
listed the difference with the analysis of Venn (1995b). 



HR 1865 

T eff =7500 K 

log 9=2.0 
it=3.0 kms" 1 
[Fe/HJ=0.00 



ion 


N 


Wx 


e 


[el/Fe] 


a 




AVenn 


C I 


11 


29 


8.16 


-0.40 





09 


+0.02 


O I 


5 


34 


8.44 


-0.49 





02 


-0.24 


Na I 


3 


67 


6.72 


+0.39 





11 


-0.20 


Mgl 


2 


23 


7.52 


-0.06 





05 


-0.26 


Si I 


9 


19 


7.34 


-0.21 





26 


-0.08 


S I 


3 


29 


7.22 


+0.01 





06 


-0.22 


Ca I 


9 


9 


6.44 


+0.08 





13 


-0.06 


Sc II 


1 


14 


3.27 


+0.17 






+0.23 


Ti I 


4 


16 


4.97 


-0.02 





07 




Ti II 


2 


101 


5.07 


+0.08 





09 


+0.19 


Cr I 


2 


15 


5.57 


-0.10 





05 




Cr II 


4 


90 


5.88 


+0.21 





12 


+0.21 


Mn I 


3 


17 


5.46 


+0.07 





13 




Fe I 


59 


36 


7.50 


-0.01 





12 


-0.15 


Fe II 


18 


50 


7.53 


+0.02 





15 


+0.10 


Ni I 


28 


17 


6.18 


-0.07 





14 


-0.08 


Cu I 


1 


5 


3.98 


-0.23 








Zn I 


1 


17 


4.36 


-0.24 








Y II 


8 


47 


2.41 


+0.17 





16 




Zr II 


5 


18 


2.98 


+0.38 





21 


+0.19 


Ce II 


5 


7 


1.71 


+0.16 





09 




Nd II 


7 


9.2 


1.73 


+0.23 
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4.1. HR 1865 and HR 1017 

Although our primary goal was to use the spectra of the 
narrow-lined supergiants HR 1865 and HR 1017 for line- 
identification purposes and to recognize possible blending 
in the spectra of the IRAS objects, we also performed 
a complete LTE abundance analysis of these two bright 
objects. 
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Table 3. Chemical analysis of HR 1017. Table 4. Chemical analysis of IRAS 22223+4327. 



HR 1017 

T eff =6500 K 
log g =2.0 
^ = 3.50 kms^ 1 
[Fe/H]=-0.16 



ion 


XT 

IN 


T T7" 

W\ 


e 


[el/rej 


a 


C I 


10 


47 


8.43 


+0.03 


0.11 


O I 


2 


32 


8.65 


-0.12 


0.14 


Na I 


1 


91 


7.01 


+0.84 




Al I 


1 


25 


6.31 


+0.00 




Si I 


8 


33 


7.19 


-0.20 


0.22 


Ca I 


7 


78 


6.24 


+0.04 


0.11 


Sc II 


3 


43 


3.19 


+0.25 


0.06 


Ti I 


7 


16 


4.84 


+0.01 


0.14 


Ti II 


4 


94 


4.83 


+0.00 


0.12 


Cr I 


7 


44 


5.50 


-0.01 


0.18 


Cr II 


7 


87 


5.75 


+0.24 


0.13 


Mn I 


6 


26 


5.32 


+0.09 


0.16 


Fe I 


69 


44 


7.32 


-0.03 


0.12 


Fe II 


12 


47 


7.39 


+0.04 


0.07 


Ni I 


30 


36 


6.22 


+0.13 


0.13 


Cu I 


1 


36 


3.98 


-0.07 




Zn I 


1 


38 


4.03 


-0.41 




Y II 


3 


37 


2.01 


-0.07 


0.11 


Zr II 


2 


33 


2.64 


+0.20 


0.25 


La II 


2 


20 


1.44 


+0.38 


0.04 


Ce II 


5 


17 


1.37 


-0.02 


0.06 


Nd II 


15 


23 


1.39 


+0.05 


0.14 



The non- variable HR 1865 is a very bright Galactic F- 
type supergiant which is often used as a standard star for 
Galactic and extragalactic research. The most recent ex- 
tensive chemical study is by Venn (1995a, b). She deduced 
model-parameters of T c ff = 7400 and log(g) = 1.1 and £ t 
= 4.0 which are similar to our findings. The difference be- 
tween the two studies is given in Table ||. Generally the 
agreement is rather good with the noticeable exception 
of Mg I, S I and Sc II where the difference amounts to 
0.26, 0.22 and 0.23 dex respectively. There are no Sc II 
and Mg I lines common in our analyses. The common S I 
lines (3 in total) have very similar equivalent widths, so 
the difference in results is due to slightly different atomic 
parameters and microturbulcnt velocity. 

The most recent analysis of HR 1017 we found in the 
literature is by Luck & Lambert (1985). The atmospheric 
parameters they found are somewhat cooler than ours : 
T cS = 6250, log(g)=0.90 and & = 3.0. Their iron abun- 
dance is somewhat higher, but generally the agreement 
between the abundances is good. 
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4.1.1. Intermediate and heavy elements (Na-Zr) 

Except for Na, the other elements do not seem to be en- 
hanced (see Fig. ||). This corresponds to the abundance 
analysis of Venn and Luck & Lambert. 

Both HR 1865 and HR 1017 display a strong sodium 
enrichment ([Na/H]=+0.39 and +0.68). Sodium enrich- 
ment in F-K supergiants was first reported by Luck (1977, 
1978). Since then several authors have detected an en- 
hancement of Na (see Venn 1995b and references therein). 
Boyarchuk & Lyubimkov (1983) proposed two possibilities 
for the Na-enrichment: 1) non-LTE effects in the tenuous 
and cool atmosphere, of which the corrections in the abun- 
dances are < —0.1 dex for weak lines 2) during a NeNa 
proton capture reaction ( 22 Ne(p,7) 23 Na) Na could have 
been synthesized and after deep mixing with the interior 
(first dredge-up), Na-enrichment may occur in the stellar 
atmosphere (Venn 1995b). The Na-synthesis can only take 
place when the temperature is high enough. Therefore Na- 



Leen Decin et al.: Optical spectroscopy of carbon-rich post-AGB stars 



7 



Table 5. Chemical analysis of IRAS 04296+3429. 
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enrichment is mass dependent: the more massive the star, 
the stronger Na will be enhanced. Takeda & Takada-Hidai 
(1994) calculated that HR 1865 is a less massive star than 
HR 1017 (M H r ioi7=17 Mq and M HR i 865 = 14.5 M Q ); 
the Na enrichment displays the same trend. 

Interestingly in an LTE analysis of several high- 
luminosity LMC/SMC Cepheids, Hill et al. (1995) did not 
found a general Na enhancement. Since the same tem- 
peratures and densities apply for the LMC/SMC Cepheid 
models, there is a hint that the Na enhancement of the 
Galactic Cepheids is real and not due to non-LTE effects. 
Detailed evolutionary models coupled with accurate nu- 
cleosynthetic networks are, however, needed in order to 
quantify the dependance of the nucleosynthetic yields and 
mixing ratios on the overall chemical composition, or other 
as yet unexplained difference in the LMC / SMC and Galac- 
tic supergiants. 

4.2. IRAS 22223+4327 and IRAS 04296+3429 

In Tables | and | we give an overview of the abundances 
relative to the solar value for IRAS 22223+4327 and 
IRAS 04296+3429. 

4.2.1. Population 

Both IRAS 22223+4327 and IRAS 04296+3429 are iron 
deficient, [Fe/H]=— 0.4 and —0.7 respectively. Together 
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Fig. 3. The abundances of HR 1865 and HR 1017 relative to 
the solar value. 



with the value of their galactic latitude (b=— 11.6 for 
IRAS 22223+4327 and b=-9.1 for IRAS 04296+3429 ) 
this indicates that the IRAS-sources are low-mass ob- 
jects of a relatively old population. The radial velocity de- 
duced from the spectra are Vlsr = —29 ± 1.0kms _1 and 
-62 ± 1.0 kms- 1 for respectively IRAS 22223+4327 and 
IRAS 04296+3429 which are similar to values found in 
the literature (Omont et al. 1993). So far their is thus no 
evidence for radial velocity variations and binarity of the 
programme stars. 

4.2.2. 3rd dredge-up 

Both stars display a strong enrichment of carbon ([C/Fe] 
= +0.5 and +0.9), an indication that the third dredge-up 
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The most convincing argument for mixing products of 
the helium-burning shell into the atmosphere are the large 
s-process elemental abundances. For IRAS 22223+4327 
and IRAS 04296+3429 we could derive abundances of 9 
respectively 7 s-process elements and all turned out to 
be significantly overabundant, even relative to the solar 
value. We will focus on the distribution of these elements 
in a separate section. 

For metal deficient stars the abundances of s-process- 
elements scale with Fe for — 1.5<[Fe/H]<0.0 (see Wheeler 
1989). For IRAS 22223+4327 and IRAS 04296+3429 
[s/Fe] is 1.0 and 1.4! There is thus no doubt that both 
IRAS stars are post 3rd dredge-up stars. 



4-3. Intermediate-mass elements (Al-Ni) 
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Fig. 4. The abundances of IRAS 22223+4327 and IRAS 
04296+3429 relative to the solar values. 



was effective. During this third dredge-up material out of 
the helium-burning shell is mixed into the stellar photo- 
sphere. Unfortunately, the O abundance is difficult to com- 
pute in this temperature-gravity domain. The only lines 
available are the O triplet at 6150 A, but these are heavily 
blended with a Fe I and a Si I line. We only could deter- 
mine the O abundance for the bright IRAS 22223+4327 
and we obtain a C/O ratio of 1.3, but due to the large 
uncertainty of the O abundance this ratio is not very ac- 
curate. From the atomic photosphcric lines alone the C/O 
ratio cannot be determined accurately enough to claim 
the objects to be real carbon stars (C/O > 1). The sum 
of the CNO nuclei is therefore also rather uncertain and 
compared to the expected value for an unevolved star of 
the same metallicity we find only a A^CNO of +0.05. 



For unevolved metal deficient stars a slight overabun- 
dance of the a-elements reflects the chemical history of 
our Galaxy. The solar neighbourhood displays an increase 
of the [a/Fe] ratio from 0.0 to +0.4 in the metallicity 
region from about solar to [Fc/H]=— 1.0 (Edvardsson et 
al. 1993) When we take into account the error bars in 
the a-abundances, the overabundances of Si, S and Ni in 
IRAS 04296+3429 and Si and Ti IRAS 22223+4327 do 
not seem to be abnormal. We may conclude that the 3rd 
dredge-up did not enhance significantly the a isotopes. 

Also the odd elements are interesting tracers of internal 
nucleosynthesis and structure. Theoretical models predict 
an enrichment of Al and especially the Al/Mg ratio in stars 
where Hot Bottom Burning (HBB) has taken place (Lat- 
tanzio et al. 1996). In the deepest layers of the convective 
envelope the temperatures may reach about 82 million K, 
and substantial hydrogen burning via the CNO cycle will 
take place. Further, at these high temperatures 26 Mg suf- 
fers substantial proton captures and produces 27 A1 during 
the interpulse phase. In theoretical evolutionary models, 
HBB is predicted only in intermediate mass stars, since 
the temperature of the bottom of the convective envelope 
only reaches in these models high enough values for the 
synthesis to take place. Although the high C abundance 
already indicates that HBB was not very effective in the 
IRAS objects, we carefully analysed the Al lines in our 
spectra of IRAS 22223+4327. The strongest optical lines 
at 8773 A (multiplet number 9) are unfortunately heavily 
blended with the Phillips (2,0) band of the circumstellar 
C2. We therefore based or analysis of the lines of multi- 
plet 10 (7835.3 and 7836.1 A) and multiplet 5 (6696.0 and 
6698.7 A). The small overabundance of [Al/Fe] = +0.2 is 
not significant. Also the 7 Li resonance line at 6707 A is 
not detected, so we conclude that there is no evidence for 
even moderate HBB in IRAS 22223+4327. Unfortunately 
there arc no useful Mg lines in our spectra. 
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4-4- s-proces element distribution 

Several theoretical studies of large nuclear reaction net- 
works coupled with accurate AGB evolutionary codes ex- 
ist (Malaney 1987a,b; Busso et al. 1992,1995) that enable 
us to characterize the s-process in the IRAS sources based 
on the photospheric s-process element distribution. The 
direct physical information on the efficiency of the internal 
nucleosythesis that can be deduced from the s-process dis- 
tribution of an individual object, is unfortunately limited 
since the predicted photospheric distribution is not only 
determined by the nucleosythesis itself, but also by the 
theoretically less understood dredge-up process (Busso et 
al. 1995). The parameters governing the outcome of such 
chemical evolutionary models do therefore not only con- 
sist of nucleosynthetic quantities, but also of more ad-hoc 
adopted values governing the stellar evolution like mass- 
loss, dredged-up mass, frequency of dredge-ups etc. (Busso 
et al. 1992). 

The ratio of light (Sr, Y, Zr) to heavy (Ba, La, Nd 
and Sm) s-elements gives a measure of the neutron expo- 
sure rate which is defined as J N n (t')V(t')dt' where N n is 
the neutron density, V the relative velocity of the neutrons 
and the seed nuclei. The pulsed neutron irradiation during 
the AGB evolution is parameterized by defining a mean 
neutron exposure tq defined a t d = At(— lnr) with At 
the neutron exposure rate of a particular pulse and r the 
overlap factor which is the fraction of the inter-shell ma- 
terial that remains in the neutron-exposed region (Ulrich 
1973). In most calculations, the neutron density and the 
overlap factor are adopted as constants and the efficiency 
of s-processing is indicated by one parameter To. 

For IRAS 22223+4327 we obtained the most complete 
s-process distribution. In the temperature-gravity domain 
of the moderately deficient programme stars, the Sr abun- 
dance is extremely difficult to measure due to the lack of 
weak lines. We therefore did not take this element into ac- 
count to determine the [ls/Fe] ratio but following Busso et 
al. (1995) no correction factor is needed to account for un- 
observed elements from the light s-process trio. The mean 
abundance of the light elements Y and Zr is [ls/Fe] = +1.5 
while for the heavy elements Ba, La, Nd and Sm [hs/Fc] = 
+0.9, hence the [hs/ls] = —0.6. Following Fig. 6 of Busso 
et al. (1995) the mean neutron exposure of the object can 
be estimated to be in between 0.2 and 0.25 mbarn -1 . 

IRAS 22223+4327 is located in this diagram in the lo- 
cus of the Carbon stars and the high value of [ls/Fe] indi- 
cates a high mixing ratio between dredge-up material and 
residual mass of the envelope (Busso et al. 1995). This is 
not surprising since the post-AGB star is Carbon rich. 

More surprising is, however, the low value of the mean 
neutron exposure of the object given its low metallicity. 
Indeed, several observational evidences exist that the neu- 
tron exposure increases with lower metallicity. In Fig. 1 of 
Busso et al. (1995), where they show the measured [hs/ls] 
values as a function of the iron abundance [Fe/H] of a 



sample of intrinsic and extrinsic S stars and Ba stars, a 
clear correlation is seen with an increase of 0.2 dex in 
[hs/ls] for a drop in metallicity from to —0.5. While the 
mean trend of the Ba-stars indicate that an object with 
metallicity of —0.5, should have a [hs/ls] of +0.2, this ratio 
is only -0.6 for IRAS 22223+4327! Since the metallicity 
of Carbon stars is difficult to measure, no observational 
material is present for comparison. 

In order to characterize more precisely the s-process 
efficiency, we made use of the abundance tables for s- 
process nucleosynthesis of Malaney (1987b). The tables 
list element abundance enhancements for exponential dis- 
tribution of neutron exposures for different values To using 
fixed neutron densities of 10 s and 10 12 cm~ 3 . To quantify 
the comparison between the observations and models we 
used the goodness-of-fit procedure as defined by Cowley 
& Downs (1980) and often used in the literature. For a 
description of the method we refer to Smith et al. (1997). 
In the comparison, we used the abundances of all the mea- 
sured species : Y (Z=39), Zr(40), Ba(56), La(57), Ce(58), 
Pr(59), Nd(60), Sm(62) and Eu(63). For the uncertainty 
on the measured abundances we added quadratically a 
fixed value of 0.3 to the internal accuracy given in Table 5, 
this to account for uncertainties on the model atmosphere 
and uncertain log gf values. Moreover, for species with 
less than 5 useful lines measured, we adopt an internal ac- 
curacy of 0.2 dex. We corrected the observed number den- 
sities for the initial abundance by adopting [s-process/Fe] 
= and a metallicity of -0.5, but since the overabundance 
is large, this correction only marginally influences the re- 
sult. For Eu, which is primarily an r-process element with 
a different chemical history than the s-process elements, 
we used the results of Woolf et al. (1995) to estimate the 
initial abundance. In this method, the quality of the fit is 
given by the quantity S 2 , with a lower value expressing 
a better fit. In Fig. || we plot our corrected abundance 
distribution together with Malaney's (1987b) model pre- 
dictions for an exponential neutron exposure with To of 
0.2 mbarn- 1 (N n = 10 8 cm" 3 ) with an S 2 value of 0.6. 
For higher neutron efficiencies, the S 2 parameter increases 
fast (e.g. for t = 0.3, S 2 is already more than tripled) and 
for the To = 0.1 model we have S 2 = 3.7. Also the model 
with t = 0.05 gives a good fit but the absolute enhance- 
ment of the s-process elements is, in this model, predicted 
to be too low in comparison with the observed values. We 
can conclude that the distribution of s-process elements 
in IRAS 22223+4327 points to a low s-process efficiency 
characterized by a To < 0.2 mbarn -1 . 

The s-process distribution of IRAS 04296+3429 also 
points to a low s-process efficiency with the best fit again 
obtained for the models with To < 0.2 mbarn -1 . In this 
object, only 7 elements were taken into account in the 
fitting. 
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Table 6. Summary of spectral features in the best studied 21 (im sources, emphasizing correlation with carbon. 
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5. Discussion 

The chemical composition indicates IRAS 22223+4327 
and IRAS 04296+3429 to be 'bona-fide' post-AGB stars: 
they display an enrichment of carbon with C/O higher 
than solar and especially high overabundances of s-process 
elements. For both stars the predicted heavy element 
abundances for an exponential distribution of neutron 
exposures with low value of the mean neutron exposure 
to = 0.2 mbarn -1 fit the observed s-process abundances 
best. 

Note that not all stars, thought to be in-between the 
AGB and PN stage of stellar evolution, have all these 
chemical properties. E.g. HD 133656 ([Fe/H]=— 1) shows 
CNO enrichment, but no overabundance of the s-process- 
elements (Van Winckel et al. 1996b). Also SAO 239853 

([Fe/H]=-0.8 1.0) (Van Winckel 1997) and HR 4912 

(Luck & Bond 1989) are thought to be post-AGB stars 
without an enrichment of the s-process-elements. On 
the contrary, other 21 yum stars like IRAS 07134+1005 
([Fe/H])=-1.17; Klochkova 1995) and IRAS 05341+0852 
([Fe/H] = -0.9; Reddy et al. 1997) do also show over- 
abundances of s-process-elements. 
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Fig. 5. Comparison between the observed s-proces abundance 
distribution of IRAS 22223+4327 and model predictions from 
Malaney (1987b) for an exponential neutron exposure char- 
acterized with ro = 0.2mbarn _1 and a neutron density of 
N n = 10 8 cm" 3 . 
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The two IRAS sources IRAS 22223+4327 and 
IRAS 04296+3429 belong to the small group of stars dis- 
playing the unidentified 21 /im emission feature. The other 
best studied 21 /im sources are listed in Table || (Kwok et 
al. 1989; Hrivnak & Kwok 1991a,b; Kwok 1993; Hrivnak 
1995; Justtanont et al. 1996). Some of the sources are re- 
solved in the 8-13 /im window (Meixner et al. 1997) Note 
that the table does not include the new possible candi- 
dates given by Henning et al. (1996). 

Table ^ shows that the 21 micron sources have some 
common properties: 

1. they are post-AGB stars with a central star of spectral 
type F or G and luminosity class I (Kwok et al. 1995) 

2. all have double-peaked energy distributions (Hrivnak 
1995) 

3. all are carbon-rich objects based on their optical, in- 
frared or millimetre spectra 

4. all seven of which have been observed show HCN molec- 
ular lines 

5. five sources display the 3.3 /im emission band and three 
of them also an unusually strong 3.4 /im emission band 

6. four display the 30 /im emission band (Hrivnak 1995; 
Omont et al. 1995) 

7. all the 21 micron sources of which the opti- 
cal spectra have been analysed (IRAS 22223+4327, 
IRAS 04296+3429, IRAS 07134+1005, IRAS 05341+0852 
and HD 187885) show an overabundance of s-process- 
elements, but also the reverse is true : all field post-AGB 
stars known to be s-process enriched, also display the 21 
/im feature. 

The carbon-rich chemistry, as evidenced by the C2 cir- 
cumstellar molecular lines and the photospheric analyses 
corroborates the suggestion of Hrivnak et al. (1989a) that 
carbon is a major constituent of the molecule produc- 
ing the 21 /im emission feature. So far the only star that 
evolved hot enough to display photospheric helium lines 
is HD 187885 (Van Winckel et al. 1996a) which only has 
a very weak 21 /im feature. The effective temperature of 
8000 K marks probably the region where the carriers of 
the feature are destroyed. High-resolution spectroscopy of 
more central stars and certainly accurate determinations 
of the effectieve temperatures are needed to test several 
theories concerning the nature of the IR features and the 
influence of changing properties of the stellar radiation 
field on them (Buss et al. 1990; Begemann et al. 1997). 

So far, the study of heavy element nucleosynthesis 
to constrain 3rd drege-up models, concentrated on two 
groups of stars : the AGB stars (mainly Carbon stars, but 
also M stars that display s-process enhancements) and the 
Ba stargroup (Ba giants, Tc poor S stars, CH stars etc.). 
In the latter group the s-process enhancement is generally 
accepted to originate from mass-transfer episodes. These 
objects are binaries with one component being a cool old 
white dwarf (WD). The chemical peculiarities were build 
up when the star, which is now the primary, accreted s- 
process element enhanced material from the companion, 



which then was an AGB star and now is a cool white 
dwarf. The s-process distribution seen in these stars is 
thus not a direct result of internal mixing of the Ba-star 
itself, but a result of a not well understood accretion pro- 
cess of s-process enhanced material. 

The study of the s-process element distribution of AGB 
stars themselves is made extremely difficult by the large 
photospheric molecular opacity that makes the modelling 
difficult and the detection of several species impossible 
(e.g. Utsumi 1985). 

A systematic study of post-AGB stars in general and 
of the 21 /im post-AGB stars in particular can in this de- 
bate contribute invaluable information. Indeed, not only 
the temperature-gravity domain of these stars makes it 
possible to detect a large number of s-elements by their 
atomic lines, but also the photospheric molecular opacity 
is negligible. Since the photospheric abundance is a direct 
result of the dredge-up process (contrary to the case of 
the Ba-stars) and a fair number of 21 /im stars are now 
known, a systematic study of the s-process distributions 
can be used directly to constrain evolutionary models. In 
this first study the s-process distribution is characerized 
for IRAS 22223+4327 and IRAS 04296+3429 both with 
mean neutron exposure r < 0.2 mbarn -1 and for an- 
other objest IRAS 05341+0852, Reddy et al. (1997) find 
the To ~ 0.3 mbarn -1 . This value is certainly much lower 
than what might be expected from their metallicity if the 
21 /im stars would follow the same metallicity-neutron 
exposure relation of the Ba-stars (see Fig. 1 Busso et al. 
1995). It is clear that a more homogeneous study of more 
sources is needed to compare carefully the different groups 
of stars. 
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